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ABSTRACT 
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upon  laser-beam  propagation  within  the  atmosphere  has  been  investigated  and 
has  been  applied  to  an  analysis  of  laser  safety. 
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STATISTICAL  MODEL  FOR  LASER  SAFETY  ANALYSIS 


1 .  INTRODUCTION 
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The  objective  of  this  study  is  to  investigate  a  statistically-based 
model  of  low-energy  pulsed  laser  propagation  within  the  atmosphere,  and  to 
apply  it  to  the  determination  of  ocular  hazard  distances. 


Safety  calculations  based  upon  the  Australian  Laser  Safety  Standard 
[1]  AS  2211  (1981)  are  primarily  geometrical  and  deterministic.  They  attempt 
to  allow  for  the  effects  of  scintillation  within  the  beam  by  applying  a 
constant  correction  factor  to  the  calculated  mean  irradiance.  However, 
ocular  hazard  distances  determined  by  this  method  may  be  conservatively  high, 
and  it  is  therefore  desirable  to  compare  the  result  with  that  obtained  using  a 
satisfactory  model  of  the  turbulence  effects. 


Procedures  for  protecting  observers  exposed  to  laser  radiation 
within  the  100  nm  to  1  mm  wavelength  range  are  described  in  AS  2211.  The 
laser  hazard  is  defined  in  terms  of  the  maximum  permissible  exposure  (MPr); 
this  is  the  quantity  of  radiant  energy  per  unit  area  specified  as  being  sate 
to  the  eye.  The  MPE  for  nanosecond-pulsed  visible  radiation  given  in  AS  2711 
is  5  x  10-3  J/m ■*,  while  for  near  infra-red  radiation  the  corresponding  value 
is  an  order  of  magnitude  higher.  These  thresholds  appear  to  have  been 
derived  from  the  1978  paper  by  Sliney  (2). 


A  nominal  ocular  hazard  distance  (UOHD)  is  described  in  the  AS  2711 
as  the  distance  at  which  the  beam  irradiance  corresponds  to  the  mpu  lor  a 
given  pulse  duration.  Using  the  'top-hat'  model  of  the  laser  bean  irradianoe 
the  expression  for  NOHD  is  normally  written  as 


NOHD 


i  /W 

$  V  "*MPE 
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where  E  is  the  laser  pulse  energy,  $  is  the  half-angle  beam  divergence,  i'; 
the  beam  homogeneity  factor  (>  1),  and  hj  is  the  lenticula tion  safety 
factor  (>  1),  The  value  =  10  is  recommended  in  AS  7211  as  a  suitable 
safety  margin  for  calculating  NOHO  when  atmospheric  turbulence  is  or  . 
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The  model  developed  in  this  report  is  statistical,  in  line  with 
previous  US  and  UK  work  involving  laser  scintillations  applied  to  safety 
clearance  determinations  [2-6J.  An  interactive  computer  program  has  been 
developed  so  that  a  range  of  input  parameter  options  can  be  made  available  for 
investigating  the  physical  processes  attending  beam  propagation. 


2.  THE  PROPAGATION  MODEL 


It  is  normally  assumed  in  laser  propagation  studies  that  atmospheric 
attenuation,  attributable  to  Rayleigh  and  Mie  scattering  and  to  water- vapour 
absorption,  follows  the  exponential  Bouguer  law  [7].  The  extinction 
coefficient  for  a  particular  wavelength  can  be  determined  using  Koschmieder's 
formula  {71  for  the  0.55  pm  wavelength: 


a0.55 


3.91 2/V 


(2) 


where  V  is  the  prescribed  visibility  {in  km).  Extinction  coefficients  at 
other  wavelengths  X  may  be  obtained  by  using  eqn  (2)  in  conjunction  with  an 
interpolation  procedure  applied  to  empirical  data.  Turbulent  refractive 
effects  associated  with  atmospheric  propagation  are  often  expressed  by  means 
of  the  refractive  index  structure  constant  Cn2  [8,9]  which  exhibits  a  height- 
dependence  often  expressed  in  the  form 
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Cn2(z)  *  Cn2(z0).lz/z0) 


(3) 


where  z  is  the  height  above  ground  level,  Zg  is  a  prescribed  height  (usually 
1  n),  and  fS  is  an  exponent  determined  by  the  environmental  conditions  (under 
the  conditions  considered  in  this  study,  it  is  valid  to  set  (J  *  1.3  18,91). 


The  hazard  ranqe  D  (in  metres)  for  a  laser  aimed  at  a  target  can  be 
estimated  from  a  deterministic  model  which  assumes  a  gaussian  irradiance 
profile  at  all  propagation  distances  [10].  This  profile  is  characteristic  of 
lasers  operating  in  the  lowest  transverse  mode  (TEMgg) .  The  maximum  (central) 
irradiance  of  a  laser  beam  is  usually  considered  in  safety  calculations,  since 
exposed  viewers  are  taken  to  be  located  in  the  most  hazardous  position  within 
the  beam.  An  effective  radius  must  be  used  even  for  irregular  beams,  hence 
the  adopted  "equivalent"  gaussian  shape  provides  a  beam  radius  where  the 
irradiance  is  1/e  of  the  maximum  value.  A  threshold  energy  T  (joules), 
directly  related  to  the  MPE  through  the  relation  T  =  MPE  x  (eye-pupil  area), 
can  be  shown  (see  the  Appendix)  to  be  given  by 

T  =  ^1*2  E  exP(~a^D)  ~  exp(-  m2w2/r02))  (4) 
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(a)  rD  -  ((wQ  +  D*)2  +  d202)1/2  (5) 

is  the  laser  spot  radius  at  the  hazard  distance  D  (metres), 

(b)  wQ  is  the  exit  beam  radius  (metres), 

(c)  +  is  the  geometric  half-angle  spread  (radians), 

(d)  0  *  X/(2xWg)  is  the  half-angle  diffraction  spread  (radians), 

(e)  m  is  the  magnification  (>  1)  associated  with  any  extra-ocular  optics, 

(f)  w  is  an  aperture  parameter  (metres)  for  the  viewing  optics  as  defined  in 
the  Appendix. 

In  the  visible  and  near-infrared  region,  diffraction  spread  is 
sometimes  an  order  of  magnitude  smaller  than  the  spread  attributable  to  other 
effects.  On  the  other  hand,  with  solid-state  lasers  inhomogeneities  in  the 
lasing  medium  may  cause  fluctuations  in  the  output  irradiance  profile,  leading 
to  larger  diffractive  effects.  where  it  can  be  established  that  the  laser  is 
operating  in  the  lowest  spatial  mode,  it  is  normally  appropriate  to  set  k. ^  *  1 
when  using  the  present  model  to  predict  the  hazard  distance.  Equation  (4) 
now  can  be  solved  easily  for  D,  using  iterative  root-finding  techniques. 

Statistical  fluctuations  in  the  irradiance  (often  called 
scintillations  or  lenticulation)  occur  in  response  to  thermal  fluctuations  in 
the  air,  through  the  mechanism  of  temperature-dependent  variations  in  the 
atmosphecic  refractive  index.  It  has  been  known  for  two  decades  [251  that 
the  atmospheric  scintillation  of  laser  beams  poses  an  additional  safety 
problem  to  that  revealed  in  a  geometrical/deterministic  safety  analysis. 
Atmospheric  turbulence  near  the  ground  can  generate  "turbulons"  (small 
turbulent  eddies)  and  associated  scintillation  spots,  or  "hot  spots",  which 
can  move  in  an  uncertain  way  within  the  laser  beam  [5,11,12,13,14],  They  can 
travel  across  the  beam  quite  readily  when  a  cross-wind  is  present,  and  are 
present  to  the  greatest  extent  when  the  vertical  temperature  gradient  is 
large.  Clear  sunshine,  a  cloudless  sky,  and  hot  ground  are  all  conducive  to 
turbulent  effects  of  this  sort.  Hue  scintillation  fluctuations  possess  a 
temporal  frequency  spectrum  ranging  as  high  as  several  hundred  hertz,  and  tend 
to  increase  in  size  with  propagation  distance  until  a  range  is  reached  at 
which  "saturation"  (levelling  out  of  the  intensity  variance)  occurs. 


One  of  the  difficulties  in  a  laser-safety  analysis  which  allows  for 
atmospheric  effects  lies  in  the  existence  of  a  non-zero  risk  of  subjecting  an 
observer  to  a  localised  irradiance  far  larger  than  would  occur  in  a  non- 
turbulent  environment  owing  to  the  possible  presence  oE  hot  spots.  The 
calculation  of  such  a  risk  is  complicated  by  the  necessity  to  allow  for 
turbulent  beam  spreading,  scintillation-saturation,  atmospheric  attenuation, 
aperture  averaging,  wind  speed,  and  other  effects.  It  is  also  necessary  to 
model  the  fluctuation  statistics  over  a  wide  range  of  turbulence 
conditions.  The  present  model  assumes  that  the  irradiance  fluctuations  obey 
modified  log-normal  statistics  [15],  although  the  method  of  calculation  can 
accommodate  a  range  of  probability  distributions  of  the  irradiance  if 
desired.  A  log-normal  distribution  arises  when  the  receiver  is  located  in  a 
turbulent  medium  and  the  scattering  cone  is  narrow,  so  that  the  received  field 


3 


nrww-v^'.VAiw.spTi  r  t’V’.-i  y 


A’.  <-VA  aaraaas  t.j;  w.  v  .  n?.?  'av* 

f'v’vV' 

O*  r 


is  the  result  of  multiplicative  effects  instead  of  additive  effects.  Detailed 
analysis  of  this  situation  is  given  by  Strohbehn  et  al  [15],  and  references 
therein.  The  irradiance  statistics  used  in  the  model  studied  in  this  report 
differs  from  a  simple  log-normal  distribution  by  the  inclusion  of  saturation 
effects  in  the  log-intensity  variance  [16-18,19],  as  well  as  aperture 
averaging  [12]. 

In  place  of  the  deterministic  eqn.  (4),  we  now  have  to  use  the 
statistical  counterpart,  which  differs  by  virtue  of  two  modifications: 


(a)  The  lenticulation  safety  factor  k 2  mus£  be  replaced  with  a  variable 
scintillation  factor  F(og),  where  the  variable  ag  is  a  modified  log-intensity 
variance.  This  variance  is  ultimately  derived  From  the  Tatarski  variance  for 
fluctuations  of  a  spherical  wave  propagating  in  the  atmosphere,  o *  [8,9]  and 
is  given  by  the  integral 


_  .7/6  11/6 

2.24  k  '  D 
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(u)  u  ( 1  — u)  du 


(6) 


where  k  *  2*/A  and  u  =  x/D,  x  being  the  distance  variable.  The  scintillation 

factor  is  defined  to  be  the  ratio  of  the  threshold  intensity  1^  to  the  mean 
intensity  <I>.  Application  of  the  theory  of  the  log-normal  distribution 
[8,15,20]  in  conjunction  with  this  definition  then  gives: 
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where 
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a  »  (in  [1  -  (2  exp(o  )  -  1)/(1  +  kw  /D)]) 
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represents  the  aperture-averaged  [12]  deviation.  Saturation  of  the  variance 
(which  occurs  when  the  turbulence  significantly  diminishes  the  beam  quality) 
has  been  observed  experimentally  since  the  late  1960's  [16,17,19-21].  An 
early  analysis  by  de  Wolf  [18]  gives  the  saturated  variance  as 


In  [2 


exp(-  0T  )], 


which  shews  reasonable  agreement  with  experimental  data  available  at  that 
time.  Another  more  recent  approach  has  been  usefully  represented  by  the 
equation  [19] 


oT  /(’ 
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the  constants  p  and  y  being  determined  (for  a  prescribed  wavelength)  by 
empirically  fitting  relevant  experimental  data  [19]  to  the  assumed  expression. 


The  saturation  effect  can  be  interpreted  in  terms  of  the  relative 
effect  of  atmospheric  turbulence  upon  the  phase  and  amplitude  of  the  beam's 
wavefront.  At  low  C^2  or  short  range  the  effect  is  almost  exclusively  a 
phase  aberration,  while  at  higher  Cn2  or  longer  ranges  the  amplitude 
fluctuations  play  the  dominant  role  and  the  beam  tends  to  spread  out  and 
ultimately  to  break  up  into  a  multitude  of  small  regions.  The  nominal  range 
at  which  the  division  between  weak  and  strong  turbulence  ocaurs  has  been  given 
[8,9]  as 
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In  eqn.  (7)  the  parameter  b  appears;  this  is  the  upper  limit  of  the 
integral  representing  the  probability  of  exceeding  the  prescribed 
threshold,  Q: 

Q  -  1  -  (2ir)“1/2  /  exp  (-  t2)dt.  (10) 

—00 

Equations  (7)  and  (10)  both  emerge  in  the  theory  of  the  log-normal 
distribution.  Series  expansions  for  the  inverse  probability  function  b(Q) , 
with  good  convergenoe  properties,  are  also  available  [22].  Note  that  large 
(small)  values  of  q(b)  correspond  with  values  for  the  scintillation  factor 
which  .-ire  less  than  (greater  than)  unity,  This  tautology  follows  directly 
from  egns.  (7)  and  (10),  and  may  be  expressed  in  the  following  way: 
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(b)  In  place  of  the 
turbulence-broadened  spot 


expression  for  r^,  eqn. 
size 


(5),  we  substitute  the 
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where  0  is  the  short-term  beamspread  to  the  n-1  point  (radians),  and 
represents  the  radial  enhancement  due  to  turbulence  of  the  observed  beam  in 
the  short-time  exposure  limit.  Note  that  rD  '  -  >,  where  rg^  i.s  the 

short-exposure  radius  defined  in  section  5(c).  The  short-term  beamspread  has 
been  evaluated  as  [23] 
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where  rQ  is  sometimes  referred  to  as  the  Fried  coherence  length,  or  lateral 
coherence  length,  of  a  spherical  wave.  It  is  given  by  [23) 


2.10  (1.455  k20  /  C  2(u)  <1-u)5/3  du)_3/5. 


The  long-term  turbulent  beam  spread  to  the  e-'  point  (radians)  is  the  radial 
enhancement  due  to  turbulence  in  the  long-time  exposure  limit,  and  is  given  by 
123) 


6it  -  0.604  ’*0/r0*  (14) 


Equations  (11)  and  (12)  for  the  infinite  gaussian  beam  have  been  modified  by 
Breaux  [23,24]  through  the  incorporation  of  new  parameters  which  serve  to 
match  empirically  the  spot  radius  calculated  via  Fresnel  integrals.  Although 
this  procedure  should  improve  the  accuracy  of  the  model,  it  has  not  been 
employed  in  the  present  calculations. 

The  equations  given  in  this  section  have  been  utilized  in  a  computer 
program  which  calculates  the  hazard  distance  for  a  given  set  of  input 
parameters.  Ttie  program  also  can  selectively  treat  any  of  the  input 
parameters  as  a  variable.  Important  input  parameters  are  listed  in  Table  1, 
and  numerical  calculations  of  hazard  distances  using  this  program  are 
discussed  in  Section  3. 

A  simplified  flow  chart  for  the  program  is  depicted  in  Figure  1. 
Being  an  iterative  calculation,  it  proceeds  via  the  large  circular  loop,  as 
shown,  until  the  hazard  distance  has  been  calculated  to  within  a  prescribed 
accuracy.  In  order  to  initiate  the  iterative  cycle,  it  is  necessary  to 
provide  a  rough  estimate  for  the  range,  and  this  is  given  by  the  solution  of 
the  deterministic  equation  (4).  Other  quantities  which  are  calculated  by  the 
program^ include  the  scintillation  factor,  the  respective  variances 
°t  '  °s  '  al  '  sflort'*term  ar*d  long-term  beam  spreads,  the  corresponding 
radii  at  the  hazard  distance,  and  the  inverse  proljabllity  variable  b. 


6 


TABLE  1 


Summary  of  the  input  parameters 


Laser  parameters  Viewer  parameters 


Propagation  parameters 


1 .  Energy 

threshold  (T) 


I.  Energy  (E) 


2.  Divergence  ($) 

(Geometric,  Jitter) 


3.  Output  beam  radius 

(wQ) 


4.  Wavelength  (A) 

5.  Height  of  laser 

U> 


2.  Probability  of 
exceeding  threshold 
(Q) 

3.  Receiver  aperture 
statistics 


4.  Viewer  statistics 

5.  Magnification  optics 


1.  Attenuation  (a^) 
(Absorption,  Scattering) 

2.  Diffraction  (0) 


3.  Refractive  index 
effects 

(Scintillation,  Beam- 
spread,  Beam  wander) 


3.  NUMERICAL  CALCULATIONS 


Calculations  of  hazard  distance  under  various  conditions  wore 
carried  out  by  means  of  the  computer  program  described  in  section  2,  using  the 
following  quantities  (Table  2)  as  fixed  parameters: 


TABLE  2 


Program  Input  Parameters 


V 

E 

T  *= 

Q 

(J)  C3 

A 


1  .06 


23.5  km 
150  mj 
2  11J 
0.01 

/2  x  10“4  rad.  (to  p~\  half-angle) 

1  .064  lim 
9  mm 

wE  =  wp  =  3.5  mm  (refer  to  Appendix) 

0.1  km-1 

2  m 
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The  value  T  ■  2  liJ  corresponds  to  the  MPE  at  wavelength  1.064  u»  with  a  pupil 
diameter  of  7  mm.  The  program  has  been  given  the  facility  for  varying  one  of 
these  input  quantities,  the  remaining  quantities  then  assuming  the  above 
values  by  default.  The  program  determines  the  hazard  distance  at  which  a 
viewer,  looking  directly  at  the  laser,  has  a  probability  Q  of  exceeding  the 
threshold  energy  T.  It  also  calculates  the  scintillation  factor,  the 
variance  o  ,  and  the  short-term  and  long-term  beam  radii. 


The  refractive  index  structure  parameter  C „  was  allowed  to  vary 
17  —10  —2/5  n  * 

over  the  range  10  to  10  (m  '  ).  As  a  rough  guide  to  the  scaling, 

10*11  corresponds  to  very  strongly  turbulent  conditions  and  10-1^  to  very  weak 

turbulence.  Figures  2  and  3  illustrate  the  results  of  the  calculations  and 

show  hazard  distances  for  different  visibilities  and  pulse  energies 

respectively.  Hazard  distances  and  scintillation  factors  in  the  present 

statistical  model  approach  the  fixed  values  from  the  corresponding 

deterministic  model  as  cn2  becomes  smaller.  With  very  large  values  of  Cn2, 

the  hazard  distance  is  dominated  by  the  effects  of  beam  spreading  and  beam 

break-up.  These  processes  ensure  that  hazard  distances  in  a  very  strongly 

turbulent  atmosphere  fall  far  short  of  those  predicted  by  a  deterministic 

model.  The  hazard  distances  and  scinti llation  factors  usually  both  peak  in 

the  moderate  turbulence  regime. 


4.  ATMOSPHERIC  TURBULENCE 


The  turbulence  level  of  the  atmosphere  can  be  estimated  if  required 
from  a  knowledge  of  the  following  measurable  meteorological  parameters. 


(a)  Insolation:  The  intensity  of  sunlight  striking  the  ground,  or 

insolation  strength,  depends  upon  the  sun’s  elevation  angle  and  the 
attenuation  of  the  atmosphere.  Vi  empirical  formula  for  Cn2  in  terms 
of  insolation  strength  can  be  determined  [6]. 

(b)  Aerosol  Level:  There  is  some  experimental  evidence  to  show  [261  that 

Cn2  can  be  expected  to  increase  as  the  meteorological  range 
(visibility)  increases  from  8  to  22  km. 

2 

(c)  Wind  Speed:  For  all  values  of  the  insolation,  Cn  falls  gradually 

with  increasing  surface  wind  speed.  Experimental  data  is  given  in 
ref.  [61  for  wind  speed  ranges  from  1.5-6  m/s. 

(d)  Terrain  Properties:  There  is  uncertain  evidence  concerning  the 

effect  of  terrain  upon  the  turbulence  level.  Bare  earth,  sand,  rock 
etc  can  be  expected  to  be  hotter  than  grass,  giving  rise  to  more 
turbulence  [6].  Cloud  cover,  time  of  day  and  season  are  clearly 
important  considerations,  since  properties  of  the  vegetation  such  as 
the  reflectance  and  transpiration  rate  are  involved. 
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(e)  Temperature:  Refractive  turbulence  does  not  appear  to  be  strongly 

related  to  the  mean  air  temperature  except  through  the  correlation 
between  insolation  level,  temperature  lapse  rate,  and  the  mean  air 
temperature  [6].  In  other  words,  the  mean  air  temperature  alone 
should  not  provide  any  indication  of  the  turbulence  level. 

(f)  Humidity:  Scintillation  at  far-infrared  wavelengths  may  be  affected 

by  absorption  of  water  vapour,  t  the  affect  is  usually  neglected  for 
most  other  wavelengths. 

(g)  Height  above  ground:  The  refractive  index  structure  constant 

decreases  with  height,  e.g.  as  in  equation  3.  Relevant  factors 
influencing  this  variation  are  macro-meteorological  phenomena  such  as 
fronts,  variations  in  terrain,  etc.,  and  micro-meteorological 
influences  such  as  the  diffusion  time  of  heat  from  the  surface  to 
higher  altitudes.  In  the  computer  program  CQ2  is  prescribed  for  the 
standard  height  1  m,  and  corresponding  values  for  other  heights  are 
found  using  eqn  (3). 


5.  DISCUSSION 


The  type  of  physical  model  that  we  have  been  considering  is  clearly 
more  satisfactory  than  deterministic  models,  but  also  possesses  the 
limitations  listed  below.  Some  of  these  have  been  discussed  previously  by 
Sliney  [2,111  and  by  Sliney  and  WolbarBht  [13]. 

(a)  Exposure  Limit:  The  first  difficulty  concerns  the  relation 

bnl-.WMii  the  aperture  over  which  the  visible  or  near-infrared  radiation  is 
averaged  and  the  exposure  limit.  It  is  relevant  to  mention  that  irradian.ee 
fluctuations  sampled  by  a  small  aperture  [e.g.,  1  mm  diameter)  can  behave 
differently  from  those  measured  by  larger  apertures,  since  the  larger  aperture 
tends  to  average  the  fluctuations  over  its  entire  area,  thereby  reducing  the 
variance  of  the  distribution.  The  MPE,  as  specified  in  AS  2211  refers  to  an 
ocular  pupil  size  of  7  mm.  This  is  usually  considered  to  be  the  largest  (or 
worst-case)  aperture  sixe,  which  corresponds  to  evening  or  dusk  conditions. 

It  would  appear  desirable  to  permit  a  higher  exposure  limit  for  daylight 
conditions  where  the  pupil  is  constricted,  but  there  is  some  debate  on  the 
issue,  and  such  a  scheme  has  not  been  implemented.  Sliney,  in  particular, 
has  provided  arguments  and  evidence  [2,11]  which  suggest  that  the  optical  gain 
of  the  eye  (the  peak  irradiance  divided  by  the  corneal  irradiance)  is  only 
slightly  affected  by  pupil  sine  in  practice.  Other  factors  relevant  to 
exposure  limits  include  the  age  of  the  viewer,  the  duration  and  repetitiveness 
of  pulses,  and  the  laser  wavelength's). 

(b)  Beam  Profile:  A  further  problem  arises  by  reason  of  the 

assumed  'gaussian*  irradiance  cross-sectional  profile  over  all  propagation 
distances.  This  is  a  convenient  shape  to  employ  since  many  lasers  operate  in 
the  lowest  transverse  mode,  however  it  must  be  recognised  that  in  highly 
turbulent  situations,  or  at  long  distances  (several  km),  the  laser  beam  is 
likely  to  consist  of  a  large  number  of  unconnected  spots  whose  Ustribution  is 
asymmetric  in  the  short  term.  At  shorter  distances  the  gaussian  profile  is 
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clearly  preferable  to  the  simpler  'top-hat*  shape.  For  output  from  the  laser 
cavity  which  does  not  already  possess  a  gaussian  shape,  we  may  use  diffraction 
theory  in  order  to  determine  the  shape  of  an  'equivalent*  gaussian*  Given  a 
circular  top-hat  irradiance  at  the  output  aperture,  the  simplest  method  is  to  * 
set  the  top-hat  radius  equal  to  the  radius  at  the  1/e  point  on  a  gaussian  of 
equal  height.  This  procedure  attempts  to  allow  for  diffraction  in  the  far 
field,  and  ensures  that  radiant  power  is  conserved.  A  more  satisfactory 
method  would  be  to  assume  that  the  irradiance  profile  by  Fraunhofer 
diffraction  possesses  an  'Airy  pattern'  or  'sombrero*  shape  [J^ (28/£)/0] z 
where  0  =  r/w  *s  the  radial  coordinate  and  £  ■  Jew.  jj  a  gaussian  curve 

is  then  easily  fitted  to  the  central  peak,  the  ring^structure  being  ignored. 

(c)  Exposure  Time:  When  considering  laser  beam  propagation  in  a 

turbulent  medium  it  is  necessary  to  distinguish  between  short-term  and  long¬ 
term  beamspread.  With  pulsed  beams,  'long-term*  implies  averaging  over  a 
large  number  of  pulses.  Large  turbulent  eddies  tend  to  deflect  the  laser 
beam  as  a  whole  (i.e.  produce  beam  wander),  whereas  3mall  eddies  mainly  tend 
to  broaden  the  be_m.  A  photograph  of  the  beam  taken  with  an  exposure  time 
much  longer  than  the  characteristic  wander  time  would  reveal  a  broadened  spot 
with  mean-square  radius  <r2  >  +  <r2en>;  here  rst  is  the  short-exposure  radius 
and  rcen  is  the  mean  centroid  of  rgfc.  This  simple  picture  is  of  course 
incorrect  when  the  turbulence  is  strong,  since  then  the  beam  profile  is 
discontinuous,  and  moreover  the  beam  wanders  only  slightly  by  comparison  with 
beams  in  lower  turbulence.  However,  in  strongly  turbulent  media  the  long- 
exposure  photograph  appears  like  the  short-term  exposure  one,  with  roughly  the 
same  radius. 

(d)  Retinal  Irradiance:  The  irradiance  at  the  retina,  I  of  a 

beam  intercepted  by  an  unaided  eye,  has  been  derived  [14!  as  a  function  of  the 
radiance  (brightness)  of  the  source  L,  the  pupil  diameter  dp,  and  the  spectral 
transmittance  of  the  ocular  media,  t(X): 


I  -  0.27  L  x(X)d  2, 

*  P 


and  this  works  adequately  for  source  angles  greater  than  5.7  mrad, 
corresponding  to  a  retinal  image  diameter  of  100  pm.  The  laser  radiance  is 
therefore  of  importance  in  laser  safety  considerations.  As  a  laser  is 
effectively  a  point  source,  the  rays  are  normally  imaged  on  the  retina  as  a 
minimal  image.  However,  atmospheric  turbulons  can  sometimes  defocus  the  rays 
to  such  an  extent  that  the  retinal  image  is  enlarged  [2!,  thereby  (hopefully) 
offsetting  the  increased  power  entering  the  pupil  when  the  eye  encounters  a 
hot  spot.  At  most,  however,  it  has  been  found  [2,11]  that  this  effect  has 
only  a  limited  impact. 


(e)  Additional  Factors:  Other  physical  processes  which  tend  to 

suppress  the  danger  to  the  retina  from  hot  spots  in  a  highly  turbulent 
atmosphere  include  beam-spreading  and  also  the  reduction  in  coherency  and  beam 
quality  associated  with  saturation  of  the  scintillations.  Beam  spreading  is 
manifested  in  the  strong  turbulence  region  as  a  large  negative  gradient  in  the 
graphs  of  the  hazard  range  vs  Cn2  (see  figures  2  and  3).  No  such  effect 
appears  in  corresponding  graphs  of  the  scintillation  factor  vs  Cn2  however. 
Based  on  experiments  designed  to  study  the  above  effects,  Sliney  et.  al. 

[2,11;  13,14]  have  come  to  the  conclusion  that  turbulence  does  not  add  to  the 
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out-of-doors  laser  hazard  as  significantly  as  was  believed  by  previous 
researchers.  Another  complicating  factor  lies  in  the  assumption  that  laser- 
induced  retinal  injury  only  depends  upon  retinal  irradiance.  This  is  not  a 
correct  assumption  for  some  cw  and  short-pulse  laser  exposures,  where  the 
injury  threshold  irradiance  diminishes  as  the  image  size  increases. 
Scintillation  effects  are  disregarded  for  cw  lasers  since  the  exposure  becomes 
averaged  over  a  number  of  scintillations. 

(f)  probability  of  injury!  The  maximum  probability  of  ocular 
injury  attending  laser  irradiation  of  an  individual  outdoors  will  be  realised 
when: 

(i)  the  individual  is  looking  in  the  direction  of  a  pulsed  laser, 

(ii)  the  individual's  eye  is  relaxed  (focused  at  infinity), 

(iii)  the  pupil  *  .  located  at  the  average  irradiance  maximum  of  the  beam, 

(iv)  the  fovea  is  directed  at  the  propagation  axis, 

(v)  the  retina  is  more  sensitive  (absorbing)  than  average, 

(vi)  the  eye  intercepts  a  significant  hot  spot, 

(vii)  laser  jitter  and  beam  divergence  are  minimal, 

(viii)  there  is  a  large  visual  range  (weather  dependent). 


The  probability  that  all  of  these  conditions  will  occur  simultaneously  at 
distances  beyond  the  NOHD  is  extremely  small.  The  probability  of  injury  is 
enhanced,  however,  when  the  individual  views  the  laser  with  the  aid  of  a 
telescope  or  pair  of  binoculars. 

Taking  into  account  the  multiplicity  of  effects  working  to  reduce 
the  significance  of  hot  spots  in  laser  safety  estimates,  and  also  takinq 
cognisance  of  calculations  involving  scintillation  statistics  by  authors  in 
the  UK  (61,  a  preliminary  value  of  Q  =  0.01  (1%)  can  be  assumed  as  a 
convenient  statistical  boundary  for  the  probability  of  hazardous  exposure  in  a 
colinear  arrangement  of  the  optical  axes  of  laser  and  viewer.  Apart  ^vom 
this  choice,  it  is  questionable  whether  the  statistical  model  employed  in  the 
present  study  is  valid  in  the  extreme  wings  of  the  range-dependent  probability 
distribution.  In  the  sample  calculations  presented  herein  for  standard  clear 
visibility  (V  =  23.5  km),  this  value  of  Q  produces  a  maximum  value  for  the 
scintillation  factor  slightly  larger  than  4.  Even  the  excessively  cautious 
value  (in  the  judgement  of  the  author)  Q  =  0.001  has  been  found  to  produce  a 
maximum  which  is  only  between  8  and  9.  It  can  be  concluded  on  this  basis 
that  the  incorporation  of  the  features  of  scintillation  saturation  and 
turbulent  beam  spreading  (amongst  others)  in  a  statistically-based  laser 
safety  analysis  will  generate  a  safety  factor  (the  maximum  scintillation 
factor)  somewhat  lower  than  the  value  k2  =  10  recommended  in  the  AS  2211. 

The  quantity  Q,  which  is  the  probability  of  exceeding  a  specified 
threshold  in  a  colinear  arrangement  of  transmitter  and  viewer,  should  not  be 
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confused  with  the  net  probability  of  receiving  ocular  injury.  In  a  realistic 
safety  estimate,  the  value  chosen  for  Q  in  practice  will  lead  to  a  net  injury 
probability  which  is  lower  than  Q  by  several  orders  of  magnitude.  This 
reduction  in  magnitude  is  the  consequence  of  allowing  for  geometrical  exposure 
factors,  as  well  as  a  number  of  biological  and  dynamic  factors  which 
(unfortunately)  may  be  difficult  to  quantify. 


It  will  be  observed  in  figs.  4a, b  that  the  hazard  ranges  and 
scintil lation  factors  are  reduced  significantly,  as  expected,  when  Q  increases 
in  magnitude.  Also  apparent  in  figs.  4a, b  is  the  expected  result  that  the 
hazard  range  is  almost  independent  of  Q  for  very  low  and  for  very  high 
turbulence.  The  latter  case  merely  confirms  that  strong  beam  spreading  and 
super-saturation  override  simple  statistical  considerations  and  tend  to 
destroy  the  effectiveness  of  hot  spots.  Lastly,  it  should  be  emphasised  that 
probabilistic  models  are  not  new  to  laser-safety  calculations  (see,  for 
example,  the  1969  paper  by  Deitz  [25]),  although  these  earlier  models 
neglected  important  features  such  as  beam  spreading.  It  is  also  clear  that 
buried  in  even  a  deterministic  approach  to  laser  safety  are  some  aspects  of 
the  probabilistic  approach,  such  as  the  relationship  between  MPE  and  the 
statistical  distribution  of  responses  from  an  ensemble  of  eyes. 


6 .  SUMMARY 


This  analysis  of  laser  safety  for  a  pulsed  laser  beam  allows  for 
atmospheric  turbulence  effects  by  modelling  the  irradiance  fluctuations  in 
terms  of  a  known  (modified  log-normal)  statistical  distribution.  The 
formulation  of  the  model  is  based  in  the  main  upon  well-known  theoretiaal 
expressions  and  empirical  algorithms  derived  from  a  considerable  quantity  of 
experimental  data. 


An  iterative  computer  program  has  been  developed  which  can 
determine,  for  a  prescribed  maximum  permissible  exposure,  either  the  spatial 
distribution  of  the  probability  of  exceeding  the  exposure  threshold  or  the 
ocular  hazard  range  corresponding  to  a  specific  probability.  Features  of 
both  UK  and  USA  statistical  models  for  determining  hazard  range  have  been 
incorpora  ted  into  the  present  program,  which  is  both  simple  to  implement  and 
fast. 
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APPENDIX 


APERTURE  PARAMETER  FOR  VIEWING  OPTICS 


Equation  (4)  introduces  an  aperture  parameter  w,  which  will  now  be 
defined  for  an  optically-aided  eye.  Let  wp  be  the  eye-pupil  radius,  wE  be 

th.  exit-pupil  radius i  and  l.t  »,  -  »c/«  “»««  »o  ls  the  coll.eter-.pertur. 
radius  (all  in  metres),  and  m  is  the  optical  magnification.  The  appropriate 
expression  for  w  can  be  derived  and  interpreted  by  reference  to  the  following 
simplified  diagram  of  the  optics s 


auxiliary  magnifying  optics 


It  is  assumed  that  the  optical  axes  of  the  laser,  the  magnifying 
optics  and  the  eye  are  all  collnear.  The  irradiance  distribution  of  an 
idealised  gaussian  beam  at  range  n  is 


I(r,D) 


A1 


where  r  is  the  radial  coordinate,  and  r0  Is  the  beam  width  (to  the  1/e 
position)  at  that  range.  The  power  transmitted  at  the  objective  (or 
collector  aperture)  is  therefore 
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wR  At  the  receiver  plane  R  the  intensity  profile  IR(r)  must  satisfy 

P  =  /  I  (r)  2lir  dr,  hence  one  easily  obtains 

c  J  o  R 


2  -a.D  -m2r2/r  2 
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If  wE  >  wR,  then  the  power  transmitted  by  the  exit  aperture  is  given  by  eqn. 
(A2),  otherwise  the  transmitted  power  is 


2irr  dr 


P  e 


“V 


_2„ 

-m  w_ 


(1 
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Equations  (A2)  and  (A4)  represent  the  power  entering  the  eye  pupil  prc  iding 
that  wp  is  at  least  as  large  as  wR  in  the  former  case  or  wE  in  the  latter 
case.  Equation  (A2)  clearly  must  be  altered  when  w^  <  wR,  and  nqn.  (A4)  whbn 
wp  <  wE>  in  both  cases  the  power  entering  the  pupil1 is 
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). 
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Hence  it  has  been  established  that  the  optiaal  power  entering  the 
pupil  of  the  eye  is  generally  given  by 
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where 


w 


min.  ( w_,  w„,  wn). 
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RESULTS 


ITERATIVE  CALCULATION  OF  OPTICAL  HAZARD  RANGE 


FIGURE  1 


Flow  chart  for  the  iterative  calculation  of  optical  hazard 
distance. 
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FIGURE  2a  Ocular  hazard  distance  of  the  naked  eye  as  a  function  of  the 

turbulence  parameter  Cn2,  for  two  different  visibility  levels: 
(a)  V  -  23.5  km  (standard  clear),  (b)  V  -  8  km  (light  haze). 
The  default  values  (Table  2)  apply  to  all  other  parameters. 


REF,  INDEX  STRUCT  .  CONST  .  (AT*'9/ 


2 

FIGURE  2b  The  scintillation  factor  F(o  )  vs  Cn  ,  corresponding  to  the  two 

visibility  values  in  Fig.  2a. 
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The  effect  of  varying  the  laser  pulse  energy  E  upon  the  hazard 
distance  for  a  6  x  42  -  aided  eyes  (a)  E  =*  200  mJ,  (ta)  150  mJ, 
(c)  100  mJ.  Default  values  in  Table  2  again  apply  to  other 
parameters. 


FIGURE  3 
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FIGURE  4a  Ocular  hazard  distance  of  an  unaided  eye,  for  various  values  of 
the  probability  of  exceeding  the  damage  threshold  T  ■  2  yjj 
(a)  Q  -  0.01,  (b)  Q  -  0.1,  (c)  Q  -  0.5,  (d)  Q  -  0.9. 


FIGURE  4b 


Scintillation  factors  corresponding  to  the  ranges  depicted  in 
Fig.  4a. 


